QCD SUM RULE APPROACH FOR THE SCALAR MESONS AS FOUR-QUARK STATES 
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We study the two point-function for the scalar mesons o, k, /o(980) and cio(980) as diquak- 
antidiquark states. We also study the decays of these mesons into irn, Kir and KK. We found that 
the couplings are consistent with existing experimental data, pointing in favor of the four-quark 
structure for the light scalar mesons. 
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It is known that the identification of scalar mesons is difficult experimentally, and that the underlying structure of 
^ \ them is not well stablished theoretically, due to the complications of the nonperturbative strong interactions. Actually, 

the observed light scalar states below 1.5 GeV are too numerous to be accommodated in a single qq multiplet. 
| /" ( ' The experimental proliferation of light scalar mesons is consistent with two nonets, one below 1 GeV region and 
another one near 1.5 GeV. If the light scalars (the isoscalars ct(500), /o(980), the isodoublet k and the isovector 
ao (980)) form an SU(3) flavor nonet, in the naive quark model the flavor structure of these scalars would be: 



a = —=(uu + dd), fo = ss, 
> ! V 2 

\ a[5 = —j={uu — dd), a,Q = ud, a$ = du, 

' v2 



(N 



X 



K = sd, K = su. (1) 



Although with this model it is difficult to understand the mass degeneracy of /o(980) and an (980), and it is hard to 
■ explain why a and K are broader than / (980) and an (980), its use is not yet discarded 0, 0, IH H, El Some alternative 
models allow a mixing between the isoscalars. However, different experimental data lead to different mixing angle 

By the other hand, the scalar mesons in the 1.3 — 1.7 GeV mass region (the isoscalars /o(1370)), /n(1500), the 
Qh. isodoublet ^(1430) and the isovector an(1450)) may be easily accommodated in an SU(3) flavor nonet. Therefore, 
theory and data are now converging that QCD forces are at work but with different dynamics dominating below and 
above 1 GeV mass. Below 1 GeV the phenomena point clearly towards an S— wave attraction among two quarks and 
two anti-quarks, while above 1 GeV it is the P— wave qq that is manifested. [To| . 

Below 1 GeV the inverted structure of the four-quark dynamics in S- wave is revealed with /n(980), an (980), k and 
cr symbolically given by |ll| 



a = udud, fo = — -=(usus + dsds), 

v2 

Oq" = dsus, a[J = —=(usus — dsds), aj = usds, 

v2 

k + = udds, k° = udus, R° — usud, kT = dsud. (2) 

This is supported by a recent lattice calculation ^^|- I n this four-quark scenario for the light scalars, the mass 
degeneracy of /o(980) and an (980) is natural, and the mass hierarchy pattern of the nonet is understandable. Besides, 
it is easy to explain why a and k are broader than fo and ao- The decays a — > 7T7t, k — > Kir and fo, ao — ► KK 
are OZI superallowed without the need of any gluon exchange, while fo — > 7T7t and ao — > t/tt are OZI allowed as it 
is mediated by one gluon exchange. Since /o(980) and ao(980) are very close to the KK threshold, the Jo (980) is 
dominated by the 7T7r state and ao(980) is governed by the r/n state. Consequently, their widths are narrower than a 
and k. 

In the four-quark scenario it is also easier to understand why, in some three-body decays of charmed mesons, the 
intermediate light scalar meson accounts for the main contribution to the total decay rate. For example, about half 
of the total decay rate of D + — > 7r"7r + 7r + , D + — > A'~7r + 7r + and Df — > ir~ir + ir + comes from D + — ► air + , D + — > ktt + 
and Df — > / (980)7r + repectively 0], while the light scalar mesons are hardly seen in the semileptonic decays 
D+ -> n-TT+g+vg, D+ -> K-tt + £ + vi and D+ -> ir-ir + i + v t 
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FIG. 1: The D + — ^ k° transition where k° is described as a sd state, a) and b) nonleptonic decay; c) semileptonic decay. 

Consider, for example, the decays D + — > R°ir + — > -ftT~7r + 7r + and Z? + — » R°i + v^ — > K^TT + £ + i'g. If S° is a quark 
anti-quark state (/e = sd) the above mentioned decays can proceed through the diagrams in Fig. 1, from where we 
see that one has two possible diagrams for the nonleptonic decay (Figs, la and lb), and only one diagram for the 
semileptonic decay (Fig. lc). However, if R° is a four-quark state (R = usud), from Fig. 2 we see that there are 
four possible diagrams for the nonleptonic decay, while there is still only one diagram for the semileptonic decay. 
Therefore, in the four-quark scenario the probability to form a scalar meson from the decay of charmed mesons, is 
much bigger for the no nlep tonic decay as compared with the semileptonic decay. A trend that seems be followed by 
the experimental data |14j. 




FIG. 2: The D + — > k transition where k is described as a usud state, a), b), c) and d) nonleptonic decay; e) semileptonic 
decay. 

In this work we use the method of QCD sum rules (QCDSR) to study the two-point function and the three-point 
function, associated with the meson decay constant and hadronic coupling contants of the scalar mesons, considered 
as four-quark states. The first evaluation of the /o(980) as a four-quark state in the QCDSR formalism was done in 
[l6| for the meson- vaccum decay constant, and in 01 f° r the hadronic coupling contants (for a review see [5). We 
extend these works by considering all the scalar mesons in the nonet and by considering different currents. 

We follow ref. and consider that the lowest lying scalar mesons are S-wave bound states of a diquark-antidiquark 
pair. As suggested in ref. |20| the diquark is taken to be a spin zero colour anti-triplet, flavour anti-triplet. Therefore, 
the (q) 2 (q) 2 states make a flavour SU(3) nonet. The corresponding interpolating fields are: 

eabc£dec(UaCj 5 d b )(u d j5Cd%), 

ulC^d^iqd^CsJ), q = u,d, (3) 

where a, 6, c, ... are colour indices and C is the charge conjugation matrix. The other members of the nonet are 
easily constructed. 

The coupling of the scalar meson S, to the scalar current js, can be parametrized in terms of the meson decay 
constant fs as .16]: 

(0\j s \S) - V^fsm^ . (4) 
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In order to compute this parameter by QCDSR, we consider the two-point correlation function 

U(q)=i / d 4 x e*-*{0\T\j a {x)jU0)]\0)- 



(5) 



In the QCD side we work at leading order and consider condensates up to dimension six. We deal with the strange 
quark as a light one and consider the diagrams up to order m s . In the phenomenological side we consider the usual 
pole plus continuum contribution. Therefore, we introduce the continuum threshold parameter sq |2lj |. In the SU(2) 
limit the /o and ao states are, of course, mass degenerate, and we get the same decay constant for them. After doing 
a Borel tranform the two-point sum rules are given by: 
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which accounts for the continuum contribution. 
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FIG. 3: The scalar meson decay constants as a function of the Borel mass. Solid: ff = f ao ; dashed: /„; dots: f a 



In the numerical analysis of the sum rules, the values used for the strange quark mass and condensates are: 
m 8 = 0.13 GeV, (qq) = -(0.23) 3 GeV 3 , (ss) = Q.S(qq), (qga.Gq) = m 2 (qq) with m 2 = 0.8 GeV 2 and (g 2 G 2 )_ = 
0.5 GeV 4 . We estimate the decay constants by using the experimental values for the smalar meson masses [l3l : 
m a = 0.5 GeV, m K = 0.8 GeV, m fo = 0.98 GeV. For the continuum thresholds we use s£ = 1.0 GeV 2 , sg = 1.2 GeV 2 , 
s f " =1.5 GeV 2 . 

From Fig. 1 we see that we get a very stable result, as a function of the Borel mass, for ff and f K . In the case 
of f g th e stability is not so good, but it is still acceptable. The problem with these sum rules, as already noticed in 
ref. |l6j |. is that the continuum contribution is only smaller than the pole contibution for small values of the Borel 
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mass (M 2 < 0.6 GeV 2 for a and n and M 2 < 0.8 GeV 2 for fo and ao). However, for these values of the Borel mass the 
perturbative contribution is smaller than the four-quark condensate contribution. The pcrturbative contribution will 
get bigger than the four-quark condensate contribution only for M 2 ~ 1.5 GeV 2 for a and k, and M 2 ~ 2. GeV 2 for 
fo and ao- Therefore, the Borel windows used is a compromise between these two restrictions. Allowing a variation 
of 0.2 GeV 2 in the continuum thresholds we arrive at the following values for the decay constants: 

/ CT = (7.5 ±1.0) MeV, f K — (1.6 i 0.3) MeV, f fo = (1.1 ± 0.1) MeV. (8) 

In order to study the scalar-pseudoscalar-pseudoscalar vertex associated with the a — > tt + tt~ , k — > K + tt~ , /o(a ) — > 
K + K~ and fo — > ir + ir~ decays, we consider the three-point function 

T»M,q) = [ d*x d*y e^* e i9 - B (0|T{ig(ar)iS(y)it(0)}|Q), (9) 



where p = p' + q, S denotes the scalar meson and Pi, Pi are the two pseudoscalar mesons in the vertex, for which we 
use the axial currents: 

ill=d al ^ 5 u a , =Ua7M75<4, = s a ~f^ 5 u a , =w a 7 M 7 5 s a . (10) 

In order to evaluate the phenomenological side we insert intermediate states for Pi, P2 and S, and we use the 
definitions in Eqs. (0J and (|llfl bellow: 

(0\&\Pi(p)) =ip»F Pi . (11) 

We obtain the following relation 

T$ en (p,p', q) = —5 9W P ' P \ 2 ^ S ? ?T 9sp 1 p 2 P'n,qu + contributions of higher resonances , (12) 

(mj, - p 2 )(m z Pi - p' )(m z P2 - q 2 ) 

where the coupling constant gsp t p 2 is defined by the matrix element 

(Pi(p')P 2 (q)\S(p)) = g SPl P 2 . (13) 

Here we follow refs. 0, 0] and work at the pion pole, as suggested in I2II for the nucleon-pion coupling constant. 
This method was also applied to the nucleon-kaon-hyperon coupling |2j l24| . to the D* — D — tt coupling |25| and to 
the J/i/j — 7r cross section [2^ . It consists in neglecting the pion mass in the denominator of Eq. (|12fl and working at 
q 2 — 0. In the QCD side one singles out the leading terms in the operator product expansion of Eq.© that mach the 
1/g 2 term. Up to dimension six only the diagrams proportional to the quark condensate times m s and the four-quark 
condensate contribute. Making a single Borel tranform to both — p 2 = —p' 2 — > M 2 we get: 

m z a — vfi%. 3 
9kK+ «- ^ F f Kf T« (e-- /M2 - e-<^) = ^ m + { -s S ))-^ { qq)M 2 



9f K + K- ', J W° (e- m ^ M -e"V ) = m 2 + { rssy + (q q }(- SS }) 

m fo~ K 6 V I 



>) M 2 ( 



qq) + ~<ss) I M 2 ( 1 - e~ s " /M " ) . (14) 



W2V V 3 

As said in the introduction, the fo — > tt + tt~ decay is mediated by one gluon exchange. Therefore, the first diagram 
proportional to 1/q 2 (at leading order) is the eight dimension condensate proportional to the mixed condensate times 
the quark condensate. Therefore, working up to dimension eight in this case we get 

9f0 ^- ^ F 2 lff ° m S e- ml/M2 ~ e~ m %' M2 ) = jgg^M . (15) 
y/07r m 2 Q - m 2 v ; 12V2M 2 V ' 

The problem of doing a sin gle Borel transformation is the fact that terms associated with the pole-continuum 
transitions are not suppressed [27|. In ref. |27j it was explicitly shown that the pole-continuum transitions have a 
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FIG. 4: The QCDSR results for the hadronic coupling constants g av + v - (squares), g KK +-K- (triangles up), gf K+K- (triangles 
down) and Qf on + n - (circles). The solid, dashed, dot-dashed and dotted lines give the linear fit to the QCDSR results respectively. 

different behaviour, as a function of the Borel mass, as compared with the double pole contribution: it grows with 
M 2 . Therefore, the pole-continuum contribution can be taken into account through the introduction of a parameter A 
in the phenomenological side of the sum rules in Eqs. I|14|) . (|15|l . by making the substitution gsp 1 p 2 ~* 9sp±p-> + AM 2 

Using F v = a/2 93 MeV, F K = 160 MeV, = 137 MeV, m K = 490 MeV and the scalar decay constant given by 
the sum rules in Eq. 10 we show, in Fig. 4 the QCDSR results for the hadronic coupling constants. We see that, in 
the Borel range used for the two-point functions, the QCDSR results do have a linear behaviour as a function of the 
Borel mass. Fitting the QCDSR results by a linear form: gsp 1 p 2 + AM 2 (which is also shown in Fig. 4), the hadronic 
couplings can be obtained by extrapolating the fit to M 2 = 0. In the limits of the continuum thresholds discussed 
above we obtain: 

5(T7r+7r - = (3.1 ± 0.5) GeV, 9reK+7r - = (3.6 ± 0.3) GeV, 
9 f0 K+K- = 9a K+K- = (1-6 ± 0.1) GeV, g fon+7r - = (0.47 ± 0.05) GeV. (16) 

The decay width of S — » P1P2 is given in terms of the hadronic coupling gsP!P 2 as: 

r(5 - PiP 2 ) = ^A(m|, m Pi ,m%), (17) 

where A(m|, rrip , Trip ) = m A s + rnj^ + Wp 2 — 1rr? s rr? Pi — 2m 2 s m 2 Pi> — 2rn 2 Pi m 2 p 2 . 

The experimental total decay width is related with a particular decay mode through: 

r^^TTTr) = -T(S -► TT+TT"), for S = a or /„ , 
3 

T{k-*Kh) = -T(k -> K + n-). (18) 

Therefore, using the experimental results: r(/ tttt) = 40 - 100 MeV .1], T(cr -> tttt) = (338 ± 48) MeV and 
T(k -> Ktt) = (410 ± 58) MeV 13] in Eqs. JTZJ) and JTHJ) above we get 

g™X- = (2-6 ± 0.2) GeV, g^ +w _ = (4.5 ± 0.4) GeV, = (1.6 ± 0.8) GeV. (19) 

Comparing Eqs. (fT^f) and ((TU|) we see that, although not exactly in between the experimental error bars, the 
hadronic couplings determined from the QCDSR calculation are consistent with existing experimental data. The 
biggest discrepancy is for gj oW + w - and this can be understood since, probably in this case, a s corrections could play 
an important role. In the case of the decay /o(ao) ~~ > K + K~ , the coupling can not be experimentally measured due 
to the unavailable phase space. 
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We have presented a QCD sum rule study of the scalar mesons considered as diquark-antidiquark states. We 
have evaluated the mesons decay contants and the hadronic couplings associated with the a — ► 7r + 7r~, k — » K + n~ , 
/o(ao) - * K + K~ and /o — > decays, using two-point and thre-point functions respectivelly. We found that the 

couplings are consistent with existing experimental data. Therefore, we consider this result as one more point in favor 
of the four-quark structure for the light scalar mesons. 
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